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Abstract In the theoretical framework of Two Higgs Doublet Model (2HDM) 
plus three right-handed neutrinos we consider a universal treatment for the mass 
matrices, aside from that the active neutrinos acquire their small mass through 
the type-I seesaw mechanism. Then, as long as a matrix with four-zero texture 
is used to represent the right-handed neutrinos and Yukawa matrices, we obtain 
a unified treatment where all fermion mass matrices have four-zero texture. We 
obtain analytical and explicit expressions for the lepton flavour mixing matrix 
PMNS in terms of fermion masses and parameters associated with the 2HDM-III. 
Further, we compare these expressions of the PMNS matrix with the most up to 
date values of masses and mixing in the lepton sector, via a likelihood test. We 
find that the analytical expressions that we derived reproduce remarkably well the 
most recent experimental data of neutrino oscillations. 


Keywords: Neutrinos; Seesaw; PMNS matrix; 2HDM-III. 


1. INTRODUCTION 

Although highly successful in terms of its phenomenological predictions, the 
Standard Model (SM) of electroweak interactions seems incomplete from 
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a theoretical view. In its present form, it is unable to predict the masses of 
fermions (leptons and quarks), or explain why there are several families of 
such particles. One of the most interesting phenomena is presented by the 
neutrino mixing, a phenomenon known as neutrino oscillation. In concordance 
with the recent work focus on neutrino physics [5], neutrino mass scale, 
corresponding Dirac or Majorana kind of fermion, and the source of Charge- 
Parity (CP) violation are unsolved questions. For that, see the experimental 
results concerning KamLAND (KL) reactor neutrinos [12, 13, 19], with 
respect to the expectations from reference Huber-Muller(HM) spectra [12, 13]. 
In each of the current high-statistics short-baseline (SBL) reactor experiments 
RENO [19, 6], Double Chooz [1] and Daya Bay [3]. In general, if neutrinos are 
massive particles and their masses are non-degenerate, it is impossible to find 
a flavour basis in which the coincidence between flavour and mass eigenstates 
holds both for charged leptons and for neutrinos. Hence, the phenomenon of 
leptonic flavour mixing is naturally appear between three charged leptons and 
three massive neutrinos. If there exist irremovable phase factors in the Yukawa 
interactions, the CP violation will naturally appear both in the quark and lepton 
sector. 

In this context, the flavour and mass generation are two concepts strongly 
intertwined. To know the flavour dynamic in models beyond the SM, we 
need to understand the flavour mechanism and mass generation arising in 
the standard theory. In this theory, the Yukawa matrices are of great interest 
because the values of its elements define to the fermion masses, as well as its 
phases factors are related with the CP violation through the mixing matrix. 

Moreover, the flavour changing currents arise from the not simultaneous 
diagonalization of Higgs and Yukawa matrices. Particularly, we will study 
the flavour dynamics through Yukawa matrices in the 2HDM-III (see therein 
references related with this model in [9]), which into the processes comes with 
flavour violation through Higgs states, that is, it allows to appear the Flavour 
Changing Neutral Currents (FCNC) mediated by Higgs fields. 

Other models like the 2HDM-III allow the FCNC [4, 14]. The difference 
between these models is in the Yukawa structure and symmetries of the Higgs 
sector as well as the possible appearance of new sources of CP violation. In 
this work, the Higgs potential preserves the CP symmetry with the Hermitian 
Yukawa matrices. 2HDM-III predicts three neutral states and a pair of charged 
states: H \ 2 , and Hf-, [15]. 

In 2HDM-III, FCNC are kept under control by imposing some texture 
of Yukawa matrices that reproduce the observed fermion masses and mixing 
angles [7]. Using texture forms allows for a direct relation between the Yukawa 
matrix elements and the parameters related with the decay widths and cross 
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section without losing the terms proportional to the light fermions masses. 
Specifically, considering a four-zero texture Yukawa matrix, one obtains in a 
natural way the Cheng-Sher ansatz for couplings flavour mix, which is widely 
used in the literature, where flavoured couplings are considered proportional 
to the involved fermion masses [8,9]. 

This work is realized in the frame of 2HDM-I1I, considering a hybrid 
treatment of the neutral leptonic sector through type-I seesaw mechanism. 
Moreover, a four-zero texture ansatz for Dirac and Majorana neutrino mass 
matrices, left and right-handed neutrinos respectively. We perform a statistical 
analysis of neutrino mixing angles using the likelihood test. 
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2. THE 2HDM AND SEESAW MECHANISM 

In order to make a minimal extension of 2HDM by introducing right-handed 
neutrinos, we need to consider six neutrino fields; three left-handed neutrinos 

v L = (v eL , u fiL , u tL ) and three right-handed neutrinos N R = ( N lR ,N 2R ,N 2R ), 

where only the left-handed fields take part in the electroweak interactions. In 
context of Two Higgs Doublet Model plus massive neutrinos, 2HDM+3v, for 
Dirac leptons the Lagrangian of Yukawa interactions has the form: 

2 

L Y =Y J (Y:iA> k N R + Y, k L<5> k l R )+H.c., (1) 

k =1 

where L = [y t ,l j is the left-handed doublet of 517(2), the index represents 
the charged leptons. The <f> t — {tit’ll ') denotes the Higgs doublets with 
( \> k = ia /lK . Finally, the 7/ with j = /,;v, are the complex Yukawa 3x3 
matrices. In flavour space, the Dirac fermion mass matrix can be written as: 

M j = -^{vi Y x + u 2 Y 1\ (2) 

where v x 2 are the vacuum expectation values (vev) associated with each of 
the Higgs doublets. In addition, these matrices can be diagonalized through a 
unitary transformation U, such that: 

U jL M P) R = ) = diag{m n ,m j2 ,mj 3 } (3) 

where Y[ = (7 /V _ Y k !U' /K are the Yukawa matrices in the mass basis, which give 
us the shape of Fermion-Fermion-Higgs couplings. 
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Here we consider that active neutrinos acquire their small mass through some 
seesaw mechanism. Hence, it is possible to write out the following hybrid 
mass term which involves both Dirac and Majorana neutrinos 

L° +M = ~^M D N R - X -v_ l M l (u l J 1(N r ) c M r N r + H.c. (4) 


In the above expression M D is the Dirac neutrino mass matrix, while M L 
and M r are symmetric mass matrices because the corresponding mass terms 
are of the Majorana type. In this case the lepton number L is not conserved. 
In order to diagonalize the hybrid Lagrangian, equation (4), we can begin by 

rewriting to L D+M as follows: 


j^D+M 



(5) 


where n L = (v l ,(N r ) c ) and 

M d+m = (M L M d M t d M r ) (6) 

is a 6x6 complex symmetric matrix and can be presented in its diagonal 
form as: 


m = U T M d+m U = diag {\ l ,\ 2 } (7) 

where U is a 6 x 6 unitary matrix. In the case that neutrino mass matrices 
satisfy the following hierarchy condition M R 3> M D 3> M L , we obtain that 

eigenvalues of M d+m matrix take the form: 

\k,M r and X 2 « M L — M d M r 'M t d . (8) 

The previous expression is known as type-(I+II) seesaw mechanism, and it is 
just the effective mass matrix of three active neutrinos. 

3. FERMION MASS MATRICES 

In general, the Dirac fermion mass matrix has an arbitrary shape, while the 
right-handed neutrino mass matrix must be symmetric, since these latter are 
Majorana particles. In particular, in this work we consider that, respectively, 
the Dirac fermion and right-handed neutrino mass matrices are represented 
with a Hermitian and complex symmetric matrix with a four-zero texture 
shape. The explicit form of these matrices are the following 
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(9) 

where 0 B = arg j £> , J and 9 C = arg\C t j . From the expressions for the Dirac 
fermion mass matrix given in equations (2) and (9), we obtain that Y k [ Yukawa 
matrices also have a shape with four-zero texture, as shown below 
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2 (10) 

where tan/3 = v 2 tv x and u 2 = v 2 + v\ = (246.22 GeV) . 

Additionally, here we consider that the left-handed neutrinos acquire 
their small mass through the type-I seesaw mechanism, which is defined as: 
M vL = M D M R l M t d . So, from the mass matrices given in equation (9) the 

M vL matrix takes the following explicit form 
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( 12 ) 


The elements of diagonal phase matrix K are defined as ip A = arg { A uL } 
and if fj = arg j B t/l |. Also, the phase factors of M vL matrix must satisfy the 
conditions 2 arg{C vL }= aig\ 'B ijL } and 2 arg{B vL ]= arg{A uL }+arg{B t/L }. 
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where the X f ’s are the eigenvalues of M f matrix and O f is a real orthogonal 

matrix. Hence, the invariants of M f matrix are 1 

From the above expressions we may express the elements of M f 
matrices in terms of its mass eigenvalues. However, they are unable to give 
us information about the possible hierarchy in the mass spectrum. Therefore, 
a matrix with the four-zero texture shape allows to have a normal or inverted 
hierarchy in the fermionic masses. This latter hierarchy only is possible for the 
left-handed neutrino masses. 


Tr {M f } - A f + B f - X f i + X f2 + A /3 , 



(14) 


3.1 The mixing matrix as function of fermion masses 

After obtaining the neutrino mass matrix through the type-I seesaw mechanism, 
let this matrix diagonalize in the context of two different scenarios, which 
depend on the mass hierarchy imposed on the neutrino mass matrix: Normal 
Hierarchy (NH) and Inverted Hierarchy (IH). 

Normal hierarchy 

The NH in the eigenvalues of M f matrix is defined as \. 3 >\. 2 >\ r Hence, the 
mass matrix parameters in terms of mass eigenvalues and the (3,3) mass matrix 
entry, take the form 


B f — X fl + X f2 + A /3 - A f , 


/ ~ 'yi 


(15) 


C 


7 



(16) 


In this expressions for the left-handed neutrinos A f = \A f \ and Bf = |Bf| 
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According with the results, we have to take Xg = — Ay with j = 1,2,3 such 
that 
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\ n > A f > A /2 for Ay , = 
A /3 > A t > Ay, for Aj 2 = 
A Aj > Ay, for Ay 3 = 


Vi 


A 


l /2 


A 


/3 


(18) 


In case of the charged leptons: A ;1 = m e ,\ l2 = m , A 73 = m T ■ The NH is evident 
by defining the adimensional parameters M f = M f / X / : ■ Also, assuming this 
hierarchical ansatz, the heaviest particle is placed in the (3,3) mass matrix 
entry. Then, it is assumed that the parameter a f A f f is very close to 1, 
therefore one can define a f = 1 — 6 f , and the mass matrix takes the expression 


M, = 


|A/iA /2 

l-s. 


A,,A 


fry 2 


A, — X n + 6. 


■«/.«/ 


l-« ■" ' 1 f> f2 




where 


— 1 Ay,) and A 2 —11 Ay 2 , 


with Ay = Ay, / Ay 3 and X f2 = lAyJ / A /3 . 


(19) 


( 20 ) 


Inverted hierarchy 

For an inverted hierarchy (IH), the relation between the eigenvalues is 
Ay-, > Ay, > Ay 3 . Analogous to NH, the mass matrix parameters are expressed 
in terms of eigenvalues as 


B f — Ay, + X f2 + Ay 3 — A f , 


( 21 ) 
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A/] A/ 2 A/3 







A, 



( 22 ) 


(23) 


According with the results, we have to take A /; = — A /; with j = 1,2,3 such 
that 


A/2 > A f° r ^/1 — |A/i|, 

A/1 ^ A A _/hr Ay2 = |Ay 2 |, 

A/2 A ; > Ayj for Ay 3 = — |Ay 3 |. 


(24) 


For neutrinos: A y/I =mv { , \ L2 =mv 2 , \ /U = mv 3 ; and for the charged 
leptons: A W1 = m e ,\ l2 = = m T . For this hierarchy, the mass matrix is 


M = 


A/iA /3 

1 - 6 , 


AyiAya 

1 - 6 , 


-Ayi + Ayj + 5y 


J(_ e f 

1-5 ?/I?/3 


h-sM 3 


1 — 6 , 


(25) 


where 


£/i — (l — +A/i), and £/ 3 = (l — 5y + A/ 3 ), (26) 

with A /3 = A /3 / A /2 and = |a a | / A /2 • 

For a normal [inverted] hierarchy in the neutrino mass spectrum the real 
orthogonal matrix that diagonalized the fermion mass matrix with four-zero 
texture, in terms of fermion masses has the form: 
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In this matrix we have 


£n[3i — 1 — 1 + wi 
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(28) 


Now the subindex f is considering as f=u,d,u,l ■ From equations (10) and 
(27) we obtain that the elements of the Yukawa matrices in the base of the mass 
7/ obey the called Cheng and Sher relation [9] 



7 J ' 



(29) 


where k,l = 1,2,3 and (xit) are complex functions of the Yukawa matrix 
parameters and the mass matrix parameter 6= which is associated with the 


2HDM. 


3.2 The flavour mixing matrix 

The flavour mixing matrix of leptons, V PMNS arises from the lack of 
correspondence between the diagonalization of the mass matrices of the 
charged leptons and left-handed neutrinos, and this is defined as: 

Vpmns = tfU v with U vl = P„p vi (30) 

Also, the lepton mixing matrix can be written as: 

V =0 T P V, 0 (31) 

y pmns r K J 
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where P v = diag[l,e ' ,e 2 j with the phases factors <T > l *t<p B / 2—6 C and 
12 — 6 h — 6 (: . Finally, the theoretical entries of the matrix V PMNS for 
th NH [IH] are given as: 


v = - 


v = - 


v = - 

fll A 


v = - 

fl2 A 


V = - 

V 3 A 




Ai®v2[I] 








D n D w 


I *> a D m 


n \> h w6iL„, 

I D n P m 


K'”l [ 3 ] '”,2 [ .,<U,2 


D n D „ m 



' D n D „ m 




v = 


V = 


V = 








P3P2111 


' ^/2^>1[3] 


' ®/3^ ) v2[l] 




'”,'« v . [ 3 ] " ! , 2 [ 1, 6 . 5 


/J ,3«3,2, 


— + ^— (>ao-o(i -g,K' - )• 

v r " 


3, 21 


(32) 


3.3 The symmetric parameterization 

In the basis where flavour eigenstates of three charged leptons are identified 
with their mass eigenstates, the flavour eigenstates of three neutrinos can be 
written as 

K % v T ) = {v A v e2 V e3 F m1 F m2 V u3 V t1 V t2 V t3 ). (33) 

As neutrinos are Majorana particles, the nine elements of PMNS lepton 
mixing matrix can be parameterized by using three rotation angles and three 
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CP-violating phases [11]. In the so called symmetrical parametrization, the 
mixing matrix has the shape [17, 18]: 
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J 12 J 23 C c 12 l3 13 L '23 c c 12‘ 5 23 c °12 J 13 C 23 C l ' 13 l '23 


(34) 

where c-. = C 056 L and .v ;; = sinOy. In this parametrization, the relation between 
flavour mixing angles and the entries of V pMNS matrix is 

, l2 | a 

V r IZ 

■ 2/3 It r 1“ • 2/) f?2 - 2/1 Z^ 3 /Tr. 

5m W 13 = P J ,5Z/7 c/,2 =—- 1 —5Z77 v 1$ = — -(35) 

i-I p I i-I p I 

1 \ y e3\ 1 re3 


From the above expressions for the mixing angles, we can conclude that 
these are exactly the same expressions that are obtained in the Standard 
parametrization [16]. In fact, the difference between the symmetric and 
standard parametrization is explicitly manifest in the CP invariants. The 
Jarlskog invariant which is used for describing the violation in conventional 
neutrino oscillations is defined as: J cp = Jm }. 


4. NUMERICAL ANALYSIS 

In this section we make a likelihood test 'y 2 with the purpose of obtaining 
the best fit point (BFP), which allows us to get the numerical values of some 
free parameters in the \ 2 function. But before, we can take advantage of the 
last experimental data reported by Planck collaboration [2] and global fits 
of neutrino oscillations data [10], All this in order to reduce the degrees of 
freedom in the analysis. 

4.1 Neutrino mass bounds 

In the three flavour context there are six independent parameters which 
govern the behaviour of neutrino oscillations: the differences of the squared 
neutrino masses, flavour mixing angles and the Dirac CP-violating phase. The 
definition of first one is Am? = m 2 — m 2 . For a normal [inverted] hierarchy 
in the neutrino mass spectrum, we can express two of the neutrino masses in 
terms of the heaviest neutrino mass, as well as Am? parameter, as: 

an d = Am,. — A777„ hl] . (36) 

"l[3] V p] 31 l 23 l "2[ 1] \ V i2\ 32 L 21 J 
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The heavy neutrino mass must satisfy the relation ml > , and can be 

considered like the only one free parameter in the above relations, since the 
oscillation parameters Am? are experimentally determined. The values for the 
parameters A m*- at BFP±1 <t, 2ct and 3 ct reported in [10] are: 

Aot 2 2 1 (i0 _5 cF 2 ) = 7.60^;J’, 7.26-7.99, 7.11 - 8.18,[Aw 2 , (lO^eF 2 )| 

= {2.48!““, 2.35 -2.59, 2.30-2.65,2.38!"“, 2.26-2.48, 2.20-2.54. 

(37) 

In the above expressions for the parameter A m\ x the upper [lower] row 
correspond to the values for a normal [inverted] hierarchy in the mass 
spectrum. Moreover, the sum of the mass of the active neutrinos must comply 
with inequality; < 0.23, for the following actual number of active 

neutrinos N - 3.15 ± 0.23 [2]. These results are inpendent of the hierarchy of 
the neutrino mass spectrum. From equations (36) and (37) the allowed ranges 
for the neutrino masses are obtained and given in the Table 1. Also it is easy 
to conclude that for both hierarchies, there is the possibility that the lightest 
neutrino could be a massless particle. 


4.2 The likelihood test x 2 


In order to verify the viability of our hypothesis of assert that all fermion mass 
matrices have the same generic shape, namely a four-zero texture, we make a 
likelihood test x 2 i n which the estimator function is defined as: 


x 2 = E 



(38) 


Table 1: Value ranges of neutrino masses, which are obtained from equations (36) and 
(37). In addition to considering the mass constraint on heavier neutrino m 2 v >Am 2 , 

jfZJ 

and the relation ©mv.< 0.23 [2]. 


Hierarchy 

mv x (10 2 eV) 

mv 2 (10 2 eV) 

mv 3 (10 2 eV) 

Am 2 .. (eV) 

Normal 

[0,7.12] 

[0,7.18] 

[0,7.25] 

[8.72xl0 _1 ,7.18] 
[8.61xl0 _1 ,7.23] 
[8.51x10-',7.30] 

[4.98,8.69] 

[4.91,8.711 

[4.84,8.74] 

BFP 

BFP± 1 a 
BFP±2a 


[0,7.32] 

[8.40xl0 _1 ,7.37] 

[4.76,8.76] 

BFP±3a 

Inverted 

[4.87,8.19] 

[4.81.8.21] 

[4.75.8.22] 

[4.96.8.23] 

[4.89.8.24] 
[4.83,8.26] 

[0,6.58] 

[0,6.64] 

[0,6.70] 

BFP 

BFP± 1 a 
BFP±2a 


[4.69,8.23] 

[4.76,8.27] 

[0,6.76] 

BFP±3a 
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Here, the superscript th states the theoretical expressions of mixing angles 
obtained from the equations (33) and (35), while the terms with superscript exp 
states the experimental data with uncertainty a^j . The experimental data for 
mixing angles considered in this analysis are given in Table 2 [10]. 

From expressions in equations (33), (35) and (36), we can see that in general 
the x 2 function depends on five free parameters yj = j. 

But with help of the analysis performed in the previous section, the heaviest 
neutrino mass is not considered like a free parameter because its numerical 
values are determined from the experimental data. Hence, the \ 2 function has 
only four free parameters. 

Now to perform the likelihood test \ 2 , we consider that the neutrino 
masses, given in the Table 1, run into the range of 2a. The values for lepton 
masses in MeV’s are [16] 
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m e =0.5109998928 ± 0.000000011, m fL 

= 105.6583715 + 0.0000035,W/77 T (39) 

= 1776.82 + 0.16.(39) 

Then, as result of the minimizing procedure of the y 2 function, for normal 
hierarchy in neutrino masses we obtain that the values of free parameters in the 
best fit point (BFP) are the following: 


<^=—6.789x10 — 17™/, <f> 2 =2.815/-™/,^ =8.355x10 — 2, 6 V 
*=3.90x10-1,171 =5.00x10- 2 eF* 2 ,.,, =1.643x10 -9. 

As mentioned above the \ 2 function depends on four free parameters and 
three physical observables. Therefore, this function has minus one degrees 


Table 2: Experimental results of neutrino mixing angles in the ranges Ict, 2a and 3a 
[ 10 ]. 


Parameter 

BFP±la 

2 a 

3a 

sin 2 0 12 ( 10 — 1 ) 

3.26 ±0.16 

2.92-3.57 

2.78-3.75 

sin 2 9 23 (10 _ 1 )[NH] 

r zT'7+0.32 

5.67—j 24 

4.14-6.23 

3.93-6.43 

sin 2 0 23 (lO _1 ) [IH] 

r '70+0.25 
/ J -0.39 

4.35-6.21 

4.03-6.40 

sin 2 (0 n )(lO- 2 )[NH] 

2.26 ± 0.12 

2.02-5.20 

1.90-2.60 

sin 2 (0 B )(lO- 2 ) [IH] 

2.29 ±0.12 

2.05-2.52 

1.93-2.65 
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of freedom, whereby we only can obtain the BFP. However, from equation 
(40) we know the numerical values for the free parameters in the BFP. So, a 
new analysis is performance fixing the CP violation phase, since this is the 
parameter less known from the experimental point of view. But, nowadays 
there are several experiments focussed on its measurement. Then, for a normal 
hierarchy in leptonic mass spectrum, we fix the value of phases 4^ and 4> 2 , 
as well as the heaviest neutrino mass m v to the values given in equation 
(40). So, the = % 2 (S/,6^) function implies one degree of freedom. This 
last choice allows us to obtain the parameter regions at different confidential 
levels. The results related to these regions are shown in Figure 1. 


4.3 The lepton mixing angles 

Here, considering the results of the above likelihood test we study the sine of 
flavour mixing angles given by equation (35), as well as the PMNS matrix. 
In Figure 2, we show the range of theoretical values obtained at ±la as 
the experimental edge values given in Table 2. One can note that for both 
6/ vssin(Q v 13 93 ), and for the b l vssin(O n l3 23 j, results are inside the region 
of 1(7. 

As an immediate result of the above likelihood test y 2 , the flavour mixing 
matrix V pMNS is numerically computed, 

In the above section we have seen that in our theoretical framework, 
2HDM+3u, where the fermion mass matrix has a four-zero texture shape. 
We can reproduce the values of oscillation parameters in a very good 
agreement with the last experimental data. The next step in this study shall 
be to investigate the phenomenological implications of these results for 
the neutrinoless double beta decay (Ov(3|3) and the CP violation in neutrino 
oscillations in matter. 



Figure 1: For normal hierarchy. In the left graph, we show the allowed region of the 
parameters & ( and 8v. 
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Figure 2: For normal hierarchy. In the left graph, we show the allowed region of the 
parameters & 1 and 8 y to 95% C.L. 


5. CONCLUSIONS 

In the theoretical framework of Two Higgs Doublet Model type III plus 
massive neutrinos (2HDM-III+3u), we shown that can be done we outlined 
a unified treatment for the fermion mass matrices in the theory. The active 
neutrinos are considered as Majorana particles and their masses are computed 
through the type-I seesaw mechanism, where the right-handed neutrinos are 
introduced in the model as a singlet under the action of the gauge group 
of the Standard Model. In such a treatment, the mass matrices of Dirac 
and right-handed neutrinos are represented with a four-zero texture ansatz, 
which implies that the mass matrix of left-handed neutrinos have also this 
shape with four-zero texture. In fact, all Dirac fermion mass matrices are 
represented with the same generic Hermitian matrix with four-zero texture 
and a normal hierarchy in the mass spectrum. Theoretical expressions were 
derived for the elements of V pMNS matrix in function of lepton masses, two 
phases <1^ and <f>_, associated with the CP violation, and two parameters Sv 
and Su which are related with the Yukawa matrices of 2HDM-III. From the 
theoretical relations of the differences of the squared neutrino masses, and 
the experimental results reported by the Planck Collaboration and neutrino 
oscillation experiments, we obtain the allowed values for the neutrino 
masses. The parameter space exploration is done by means of likelihood 
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test x 2 ; this allowed us to find the allowed regions of the parameters 5u and 
fyat 70% and 95% C.L. for a normal hierarchy, as well as, the best fit point 
(BFP), and the mixing matrix V pMNS at 70% C.L. Finally, it is observed that 
the mixing angle as function of Sv and 8, are in very good agreement with 
experimental data. 
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